Metabolomics as a tool to explore the staphylococcal biofilm

Orthopaedic infections can be polymicrobial existing as a microbiome. Infections often incorporate staphylococcal species, including Staphylococcus aureus. Such infections can lead to life threatening illness and implant failure. Furthermore, biofilm formation on the implant surface can occur, increasing pathogenicity, exacerbating antibiotic resistance and altering antimicrobial mechanism of action. Bacteria change dramatically during the transition to a biofilm growth state: phenotypically; transcriptionally; and metabolically, highlighting the need for research into molecular mechanisms involved in biofilm formation. Metabolomics can provide a tool to analyse metabolic changes which are directly related to the expressed phenotype. Here, we aimed to provide greater understanding of orthopaedic infection caused by S. aureus and biofilm formation on the implant surface. Through metagenome analysis by employing: implant material extraction; DNA extraction; microbial enrichment; and whole genome sequencing, we present a microbiome study of the infected prosthesis to resolve the causative species of orthopaedic hip infection. Results highlight the presence of S. aureus as a primary cause of orthopaedic infection along with Enterococcus faecium and the presence of secondary pathogen Clostridium difficile. Although results were hindered by the presence of host contaminating DNA even after microbial enrichment, conclusions could be made over the potential increased pathogenicity caused by the presence of a secondary pathogen and highlight method and sample preparation considerations when undertaking such a study. Following this finding, studies were focused on an orthopaedic clinical isolate of S. aureus and a metabolome extraction method for staphylococcal biofilms was developed using cell lysis through bead beating and solvent metabolome extraction. The method was found to be reproducible when coupled with liquid chromatography-mass spectrometry (LC-MS) and bioinformatics, allowing for the detection of significant changes in metabolism between planktonic and biofilm cultures to be identified and drug mechanism of actions (MOA) to be studied. Metabolomics results highlight significant changes in a number of metabolic pathways including arginine biosynthesis and purine metabolism between the two cell populations, evidence of S. aureus responding to their changing environment, including oxygen availability and a decrease in pH. Focused investigations on purine metabolism looking for biofilm modulation effects were carried out. Modulation of the S. aureus biofilm phenotype was observed through the addition of exogenous metabolites. Inosine increased biofilm biomass while formycin B, an inosine analogue, showed a dispersal effect and a potential synergistic effect in biofilm dispersal when coupled with gentamycin. Changes in metabolism between planktonic cells and biofilms highlight the requirement for antimicrobial testing to be carried out against planktonic cells and biofilms. Untargeted metabolomics was used to study the MOA of triclosan in S. aureus. The triclosan target and MOA in bacteria has already been characterised, however, questions remain over its effects in bacteria. Although the use of triclosan has come under increasing speculation, its full effects are still largely unknown. Results show that triclosan can induce a cascade of detrimental events in the cell metabolism including significant changes in amino acid metabolism, affecting planktonic cells and biofilms. Results and conclusions provide greater understanding of orthopaedic infections and specifically focus on the S. aureus biofilm, confirming S. aureus as a primary cause of orthopaedic infection and using metabolomic analysis to look at the changing state of metabolism between the different growth states. Metabolomics is a valuable tool for biofilm and drug MOA studies, helping understand orthopaedic infection and implant failure, providing crucial insight into the biochemistry of bacteria for the potential for inferences to be gained, such as the MOA of antimicrobials and the identification of novel metabolic drug targets.

Unravelling metabolism of Leishmania by metabolomics

The leishmaniases are neglected tropical diseases with an urgent need for effective drugs. Better understanding of the metabolism of the causative parasites will hopefully lead to development of new compounds targeted at critical points of the parasite’s biochemical pathways. In my work I focused on the pentose phosphate pathway of Leishmania, specifically on transketolase, sugar utilisation, and comparison between insect and mammalian infective stages of the parasites. The pentose phosphate pathway (PPP) is the major cellular source of NADPH, an agent critical for oxidative stress defence. The PPP uses glucose, reduces the NADP+ cofactor and produces various sugar phosphates by mutual interconversions. One of the enzymes involved in this latter part is transketolase (TKT). A Leishmania mexicana cell line deleted in transketolase (Δtkt) was assessed regarding viability, sensitivity to a range of drugs, changes in metabolism, and infectivity. The Δtkt cell line had no obvious growth defect in the promastigote stage, but it was more sensitive to an oxidative stress inducing agent and most of the drugs tested. Most importantly, the Δtkt cells were not infective to mice, establishing TKT as a new potential drug target. Metabolomic analyses revealed multiple changes as a consequence of TKT deletion. Levels of the PPP intermediates upstream of TKT increased substantially, and were diverted into additional reactions. The perturbation triggered further changes in metabolism, resembling the ‘stringent metabolic response’ of amastigotes. The Δtkt cells consumed less glucose and glycolytic intermediates were decreased indicating a decrease in flux, and metabolic end products were diminished in production. The decrease in glycolysis was possibly caused by inhibition of fructose-1,6-bisphosphate aldolase by accumulation of the PPP intermediates 6-phosphogluconate and ribose 5-phosphate. The TCA cycle was fuelled by alternative carbon sources, most likely amino acids, instead of glucose. It remains unclear why deletion of TKT is lethal for amastigotes, increased sensitivity to oxidative stress or drop in mannogen levels may contribute, but no definite conclusions can be made. TKT localisation indicated interesting trends too. The WT enzyme is present in the cytosol and glycosomes, whereas a mutant version, truncated by ten amino acids, but retaining a C-terminal targeting sequence, localised solely to glycosomes. Surprisingly, cells expressing purely cytosolic or glycosomal TKT did not have different phenotypes regarding growth, oxidative stress sensitivity or any detected changes in metabolism. Hence, control of the subcellular localisation remains unclear as well as its function. However, these data are in agreement with the presumed semipermeable nature of the glycosome. Further, L. mexicana promastigote cultures were grown in media with different combinations of labelled glucose and ribose and their incorporation into metabolism was followed. Glucose was the preferred carbon source, but when not available, it could be fully replaced with ribose. I also compared metabolic profiles from splenic amastigotes, axenic amastigotes and promastigotes of L. donovani. Metabolomic analysis revealed a substantial drop in amino acids and other indications coherent with a stringent metabolic response in amastigotes. Despite some notable differences, axenic and splenic amastigotes demonstrated fairly similar results both regarding the total metabolic profile and specific metabolites of interest.

Filling in gaps of Drosophila melanogaster urate degradation metabolic pathway using metabolomics approaches: towards the core metabolome of the fruit fly

The primary goal of systems biology is to integrate complex omics data, and data obtained from traditional experimental studies in order to provide a holistic understanding of organismal function. One way of achieving this aim is to generate genome-scale metabolic models (GEMs), which contain information on all metabolites, enzyme-coding genes, and biochemical reactions in a biological system. Drosophila melanogaster GEM has not been reconstructed to date. Constraint-free genome-wide metabolic model of the fruit fly has been reconstructed in our lab, identifying gaps, where no enzyme was identified and metabolites were either only produced or consume. The main focus of the work presented in this thesis was to develop a pipeline for efficient gap filling using metabolomics approaches combined with standard reverse genetics methods, using 5-hydroxyisourate hydrolase (5-HIUH) as an example. 5-HIUH plays a role in urate degradation pathway. Inability to degrade urate can lead to inborn errors of metabolism (IEMs) in humans, including hyperuricemia. Based on sequence analysis Drosophila CG30016 gene was hypothesised to encode 5- HIUH. CG30016 knockout flies were examined to identify Malpighian tubules phenotype, and shortened lifespan might reflect kidney disorders in hyperuricemia in humans. Moreover, LC-MS analysis of mutant tubules revealed that CG30016 is involved in purine metabolism, and specifically urate degradation pathway. However, the exact role of the gene has not been identified, and the complete method for gap filling has not been developed. Nevertheless, thanks to the work presented here, we are a step closer towards the development of a gap-filling pipeline in Drosophila melanogaster GEM. Importantly, the areas that require further optimisation were identified and are the focus of future research. Moreover, LC-MS analysis confirmed that tubules rather than the whole fly were more suitable for metabolomics analysis of purine metabolism. Previously, Dow/Davies lab has generated the most complete tissue-specific transcriptomic atlas for Drosophila – FlyAtlas.org, which provides data on gene expression across multiple tissues of adult fly and larva. FlyAtlas revealed that transcripts of many genes are enriched in specific Drosophila tissues, and that it is possible to deduce the functions of individual tissues within the fly. Based on FlyAtlas data, it has become clear that the fly (like other metazoan species) must be considered as a set of tissues, each 2 with its own distinct transcriptional and functional profile. Moreover, it revealed that for about 30% of the genome, reverse genetic methods (i.e. mutation in an unknown gene followed by observation of phenotype) are only useful if specific tissues are investigated. Based on the FlyAtlas findings, we aimed to build a primary tissue-specific metabolome of the fruit fly, in order to establish whether different Drosophila tissues have different metabolomes and if they correspond to tissue-specific transcriptome of the fruit fly (FlyAtlas.org). Different fly tissues have been dissected and their metabolome elucidated using LC-MS. The results confirmed that tissue metabolomes differ significantly from each other and from the whole fly, and that some of these differences can be correlated to the tissue function. The results illustrate the need to study individual tissues as well as the whole organism. It is clear that some metabolites that play an important role in a given tissue might not be detected in the whole fly sample because their abundance is much lower in comparison to other metabolites present in all tissues, which prevent the detection of the tissue-specific compound.

The metabolic and epigenetic effects of the isocitrate dehydrogenase-1 mutation in glioma

Cancer cells have been noted to have an altered metabolic phenotype for over ninety years. In the presence of oxygen, differentiated cells predominately utilise the tricarboxylic acid (TCA) cycle and oxidative phosphorylation to efficiently produce energy and the metabolites necessary for protein and lipid synthesis. However, in hypoxia, this process is altered and cells switch to a higher rate of glycolysis and lactate production to maintain their energy and metabolic needs. In cancer cells, glycolysis is maintained at a high rate, even in the presence of oxygen; a term described as “aerobic glycolysis”. Tumour cells are rapidly dividing and have a much greater need for anabolism compared to normal differentiated cells. Rapid glucose metabolism enables faster ATP production as well as a greater redistribution of carbons to nucleotide, protein, and fatty acid synthesis, thus maximising cell growth. Recently, other metabolic changes, driven by mutations in genes related to the TCA cycle, indicate an alternative role for metabolism in cancer, the “oncometabolite”. This is where a particular metabolite builds up within the cell and contributes to the tumorigenic process. One of these genes is isocitrate dehydrogenase (IDH) IDH is an enzyme that forms part of the tricarboxylic acid (TCA) cycle and converts isocitrate to α-ketoglutarate (α-KG). It exists in three isoforms; IDH1, IDH2 and IDH3 with the former present in the cytoplasm and the latter two in the mitochondria. Point mutations have been identified in the IDH1 and IDH2 genes in glioma which result in a gain of function by converting α-KG to 2-hydroxyglutarate (2HG), an oncometabolite. 2HG acts as a competitive inhibitor of the α-KG dependent dioxygenases, a superfamily of enzymes that are involved in numerous cellular processes such as DNA and histone demethylation. It was hypothesised that the IDH1 mutation would result in other metabolic changes in the cell other than 2HG production, and could potentially identify pathways which could be targeted for therapeutic treatment. In addition, 2HG can act as a potential competitive inhibitor of α-KG dependent dioxygenases, so it was hypothesised that there would be an effect on histone methylation. This may alter gene expression and provide a mechanism for tumourogenesis and potentially identify further therapeutic targets. Metabolic analysis of clinical tumour samples identified changes associated with the IDH1 mutation, which included a reduction in α-KG and an increase in GABA, in addition to the increase in 2HG. This was replicated in several cell models, where 13C labelled metabolomics was also used to identify a possible increase in metabolic flux from glutamate to GABA, as well as from α-KG to 2HG. This may provide a mechanism whereby the cell can bypass the IDH1 mutation as GABA can be metabolised to succinate in the mitochondria by GABA transaminase via the GABA shunt. JMJ histone demethylases are a subset of the α-KG dependent dioxygenases, and are involved in removing methyl groups from histone tails. Changes in histone methylation are associated with changes in gene expression depending on the site and extent of chemical modification. To identify whether the increase in 2HG and fall in α-KG was associated with inhibition of histone demethylases a histone methylation screen was used. The IDH1 mutation was associated with an increase in methylation of H3K4, which is associated with gene activation. ChiP and RNA sequencing identified an increase in H3K4me3 at the transcription start site of the GABRB3 subunit, resulting in an increase in gene expression. The GABRB3 subunit forms part of the GABA-A receptor, a chloride channel, which on activation can reduce cell proliferation. The IDH1 mutation was associated with an increase in GABA and GABRB3 subunit of the GABA-A receptor. This raises the possibility of GABA transaminase as a potential therapeutic target. Inhibition of this enzyme could reduce GABA metabolism, potentially reducing any beneficial effect of the GABA shunt in IDH1 mutant tumours, and increasing activation of the GABA-A receptor by increasing the concentration of GABA in the brain. This in turn may reduce cell proliferation, and could be achieved by using Vigabatrin, a GABA transaminase inhibitor licensed for use in epilepsy.

Antiplatelet therapy and clinical outcomes in cardiovascular diseases

Cardiovascular diseases (CVD) is a leading cause of death in the world. Despite effective treatment regimens for ischaemic heart disease (IHD) and ischaemic stroke, mortality and recurrence rates remain high. Antiplatelet therapy is on effective treatment and reduces the risk of recurrent heart attack and stroke. Nevertheless, there are patients who stopped or interrupted their antiplatelet therapy for certain reasons or some patients may be resistant or poor responders to antiplatelet therapy. Furthermore, there is evidence of rebound effect in platelet activity after antiplatelet cessation and this may associate with increased risk of cardiovascular event. This thesis is divided into five main chapters (chapters 3 to 7) which attempt to provide data to help resolve the uncertainty. Chapter 1 highlights the background of cardiovascular diseases and the global burden of cardiovascular and cerebrovascular diseases. The metabolism of platelets, antiplatelet therapy and current antiplatelet therapy guidelines are described, followed by discussion of the risk of cardiovascular event and changes in antiplatelet therapy. Chapter 2 describes the data source from Virtual International Stroke Trial Archive (VISTA) and National Health Service Greater Glasgow and Clyde (NHSGGC) Safe Haven, followed by definition of outcome measures. In chapter 3, Virtual International Stroke Trial Archive (VISTA) data was examined to test whether continue with the same antiplatelet therapy or changing to a new antiplatelet regimen reduces the risk of subsequent events in patients who experience a stroke whilst taking antiplatelet therapy. The findings indicate that subjects who switch to a new antiplatelet regimen after stroke did not have a lower early recurrence rate than subjects who continued with the same antiplatelet therapy. Observations on bleeding complications were similar in both groups. However, changing antiplatelet regimen after stroke was associated with more favourable functional outcome across a full scale modified Rankin Scale (mRS) at 90 days. In chapter 4, association between early or later initiation of antiplatelet with a recurrent ischaemic stroke and bleeding complications was assessed using VISTA data. The findings indicate that there was no association between a recurrent ischaemic stroke and timing of initiation of antiplatelet drug after stroke. However, early initiation was associated with increased risk of bleeding. In terms of functional outcomes, this study demonstrated that the mid-time and late initiation of antiplatelet therapy after acute stroke are associated with better functional outcomes compared with early initiation. In chapter 5, a nested case-control study was performed to explore the rate of antiplatelet cessation and interruption in a sample of patients with recent ischaemic stroke and to assess the risk of cardiovascular events associated with cessation and interruption of antiplatelet. It was found that there was no increased risk of cardiovascular event among patients who had early cessation or interrupted/stopped antiplatelet therapy within 90 days following acute ischaemic stroke. In chapter 6, the incidence and predictors of cardiovascular events after DAPT cessation were evaluated. The incidence of cardiovascular event while taking DAPT and following discontinuation of DAPT was 15.7% and 16.7% respectively. This study found that increasing age was associated with an increased risk of cardiovascular event, whereas, revascularization-treated patients and longer duration of DAPT, were each associated with a decreased risk. The duration of DAPT six months and less was associated a significantly higher risk for cardiovascular event. In chapter 7, an untargeted metabolomics analysis was performed while on DAPT (aspirin plus ticagrelor) and once they stopped ticagrelor to identify metabolite changes associated with cardiovascular events after stopping DAPT. Ten ACS patients were recruited in this study and data were analysed for seven patients. Three hundred eleven putative metabolites were identified. This study found 16 putative metabolites significantly altered following ticagrelor cessation. Of these, seven metabolites were from lipid pathway and down-regulated some up to 3-fold. On the other hand, adenosine, from nucleotide metabolism was upregulated up to 2.6-fold. It concluded that there are changes in numerous pathways following DAPT discontinuation and whether these changes differ in patients who have cardiovascular event after stopping DAPT warrant further investigation. In chapter 8, a summary of the findings of this thesis are presented as well as the future directions of research in this area.